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The dynamic scaling properties of the one dimensional Burgers equation are expected to change 
with the inclusion of additional conserved degrees of freedom. We study this by means of 1-D 
driven lattice gas models that conserve both mass and momentum. The most elementary version 
of this is the Arndt-Heinzel-Rittenberg (AHR) process, which is usually presented as a two species 
diffusion process, with particles of opposite charge hopping in opposite directions and with a variable 
passing probability. From the hydrodynamics perspective this can be viewed as two coupled Burgers 
equations, with the number of positive and negative momentum quanta individually conserved. We 
determine the dynamic scaling dimension of the AHR process from the time evolution of the two- 
point correlation functions, and find numerically that the dynamic critical exponent is consistent 
with simple Kardar- Par isi- Zhang (KPZ) type scaling. We establish that this is the result of perfect 
screening of fluctuations in the stationary state. The two-point correlations decay exponentially 
in our simulations and in such a manner that in terms of quasi-particles, fluctuations fully screen 
each other at coarse grained length scales. We prove this screening rigorously using the analytic 
matrix product structure of the stationary state. The proof suggests the existence of a topological 
invariant. The process remains in the KPZ universality class but only in the sense of a factorization, 
as (KPZ)^. The two Burgers equations decouple at large length scales due to the perfect screening. 

PACS numbers: 64.60.Ht, 05.40-a, 05.70.Ln, 44.10.+i 



INTRODUCTION 



Many non-equilibrium driven systems display scale in- 
variance in their stationary states, i.e., strongly corre- 
lated collective structures without a characteristic length 
scale limiting the fluctuations. Such correlations typi- 
cally decay as power laws with critical exponents that 
are universal. Their values depend only on global issues 
like dimensionality, symmetry, and specific microscopic 
conservation laws. The classification of dynamic univer- 
sality classes and the determination of their scaling di- 
mensions is one of the central issues in current research of 
non-equilibrium statistical physics [l|, Q . The one-species 
asymmetric exclusion processes (ASEP) serves in this 
context as both the simplest prototype model for driven 
one-dimensional (ID) stochastic particle flow [1,[1,|5|] and 
as a fully discretized version of the ID Burgers equation 
(with time and space discretized, and momentum quan- 
tized) 0. 

In this paper we investigate how the properties of such 
stochastic flows change with the introduction of addi- 
tional bulk conservation laws. The generic expectation is 
that enforcing more conservation laws changes the scal- 
ing dimensions. We follow a bottom-up approach. An 
example of a top-down approach is the current interest 
in anomalous ID heat conduction in Fermi-Past a-Ulam 
type models (e.g., a chain of anharmonic oscillators 0, 
or a one dimensional gas of particles in a narrow channel 
with different types of interactions Q). The systems are 
coupled to heat reservoirs on either end. Those are held 
at different temperatures and thus induce heat flow along 
the channel. Computer simulations, e.g., using molecu- 
lar dynamics, show an anomalous thermal conductivity, 

Jq :^ /^VT, diverging with system size L diS tz ^ L^. 



The numerical estimates for the value of a in the various 
versions of the process vary between 0.22 < a < 0.44 
B Q- <^ is expected to be universal. From the ana- 
lytic side, a mode-coupling treatment predicted a = 2/5 
[lO], while a renormalization analysis of the full hydro- 
dynamic equations predicts a = 1/3, based on Galilean 
invariance and an assumption of local equilibrium in the 
heat sector [llj. In our study we add conservation laws 
to the Burgers equation instead of coarse graining down 
from full hydrodynamics. 

The equivalences between ASEP, KPZ growth, and the 
Burgers equation are well known ASEP is usually 
interpreted as a process for stochastic particle transport, 
while the Burgers equation 
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represents the evolution of a (vortex free) velocity field 
vix^t)^ and conserves momentum only [12]. The inter- 
pretation of ASEP as a fully discretized Burgers equation 
poses some conceptional issues. Due to the full quantiza- 
tion of the momenta in ASEP, in units of n = 0, 1, it can 
appear that the process also conserves energy. A careful 
discussion [13] shows that energy is conserved between 
updates but fluctuates during each update. Therefore 
ASEP is a genuine fully discretized implementation of the 
Burgers equation from this direct point of view as well. 
In section [III we discuss how to impose conservation of 
particles in addition to conservation of momentum. This 
leads naturally to the two-species ASEP known as the 
Arndt-Heinzel-Rittenberg (AHR) model [H, [TH- This 
process has been the focus of intensive studies, but its 
dynamic scaling properties seem to have been ignored. 
Instead, the stationary state properties have been cen- 
ter stage, in particular its clustering, and that it can be 
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constructed exactly using the so-called matrix product 
ansatz method [11, [13, [11, [H, Q • 

We establish that the introduction of this additional 
conservation law to ASEP does not change the univer- 
sality class, but it does so in a rather intricate manner. 
KPZ scaling changes to (KPZ)^ type scaling. The AHR 
process can be interpreted as a coupled Burgers and dif- 
fusion equation, conserving both mass and momentum; 
or as two coupled Burgers equations, one for positive and 
negative momentum quanta separately. The latter point 
of view turns out to be the most productive. At large 
length scales the coupling vanishes and the process fac- 
tor izes, in terms of quasi-particles, into two decoupled 
Burger processes. This is achieved by means of perfect 
screening of fluctuations in the stationary state. We ob- 
serve this numerically from the behavior of the two-point 
correlators (sections IIVI and |Vj) . The stationary state of 
the model is known to satisfy the so-called Matrix Prod- 
uct ansatz We use that property to prove analyti- 
cally that the perfect screening is rigorous (section rVl|) . 
In sections HVl and IVl we present also direct numerical ev- 
idence that the dynamic critical exponent is indeed the 
same as in KPZ, z = 3/2, using the time evolution of the 
two point correlators. The conventional methods fail due 
to time oscillations. This might be the first example of 
such a numerical dynamic analysis in terms of correlation 
functions. 



II. THE AHR MODEL 

Our aim is to construct a generalization of ASEP de- 
scribing a process where particle diffusion and the Burg- 
ers equation are coupled to each other. Energy will not 
be conserved. The particles in such a model need to carry 
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FIG. 1: Two species asymmetric exclusion process (bottom) 
and its corresponding interface growth model (middle) and 
particle flow model with momentum conservation (top). 



an internal degree of freedom, representing momentum. 
A site could be in four states. It would be empty = 0) 
or be occupied by a particle (ux = 1) with momentum 
Vx = +1,0,-1. Particles with +1 (-1) momentum would 
hop with a right (left) bias. Some reflection on the na- 
ture of the passing processes (the collisions) shows that 
we can remove the zero momentum state of particles, 
without loss of generality [l^ . 

This then leads naturally to the two-species ASEP 
known as the Arndt-Heinzel-Rittenberg (AHR) model. 
The conventional interpretation of this process is in terms 
of diffusion of charged particles in an electric field. Two 
species of particles with opposite unit charge hop in op- 
posite directions along a ID lattice ring, driven by the 
electric field. 
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Each site x can be in 3 states, Sx = 0,±1, with S = 1 
{S = —1) representing the right (left) moving species 
and 5 = an empty site, p is the free directed hopping 
rate (the electric field) and r the passing rate of opposite 
charged particles. In our study, the numbers of 6* = 1 
and = — 1 particles on the ring are chosen to be equal. 
Compared to the conventional single species ASEP, this 
process has two local conservation laws instead of one; 
both species are conserved independently. 

In the coupled diffusion-Burgers equation interpreta- 
tion of the same process, the charge represents a quan- 
tum of momentum moving in opposite direction as illus- 
trated in Fig ID No driving force is present, because the 
preferred hopping direction represents the total deriva- 
tive in the Navier-Stokes equation, just as in the single 
species ASEP. Similar to ASEP, energy is not a conserved 
quantity: The energy of particles is conserved between 
updates but fluctuates during the updates. That leaves 
particles in different places than where they would have 
been if energy were conserved . 

The AHR model reduces to the Sx = spin (momen- 
tum quanta) representation of ASEP in the high density 
limit where vacant sites S = are absent. There, the par- 
ticle density can not fluctuate anymore, and the process 
falls thus back to the Burgers equation with only one con- 
servation law. This limit is singular. The AHR process 
is not the generic S = 0, ±1 generalization of ASEP in 
the sense of the KPZ and Burgers equation. The proper 
generalization would be the so-called restricted solid-on- 
solid (RSOS) model (Kim-Kosterlitz model) where + and 
— pairs can be annihilated and created. Those processes 
conserve momentum. The S = {S = —1) particles 
represent up (down) steps in the KPZ type interface, the 
free hopping rate p represents step-flow. Growth at flat 
terraces is blocked in the AHR process, except for the de- 
position of vertical dimers (with rate r) in single particle 
puddles. Figil] illustrates this. 

This means that from the KPZ point of view the AHR 
process represents a growing interface where the num- 
ber of up and down steps are individually preserved. 
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Whether this local conservation law changes the scal- 
ing dimensions on large scales is the central issue we 
address here. From the KPZ perspective, your initial 
guess would probably be "no", and from the lattice gas 
perspective "yes". Our results presented below confirm 
the "no", but in rather subtle manner, the universality 
class is "(KPZ)2" instead of simple KPZ. 

The AHR model has been widely studied recently, with 
as focus the structure of its stationary state [13, EEI, HB, 
[13, El, . We are not aware of any previous dynamic 
scaling analysis. The stationary state shows strong clus- 
tering, as function of decreasing passing versus free hop- 
ping probability, r/p. Stretches of "empty" road are fol- 
lowed by high density clusters. These are mixtures of 
5* = +1 and = — 1 particles. We will identify the 
amount of mixing with the quasi-particles and the clus- 
ter size with the screening length. 

The passing of + and — particles resembles collisions. 
The ratio r/p controls the duration of the collision (the 
softness of the balls). This passing delay creates queuing 
and is the origin of the clustering. The full AHR model 
includes a reverse-passing probability (( — h) (H — ); 
particles switching position in the direction opposite to 
the electric field). That enhances the clustering even 
more. We limit ourselves here to the t = version of 
the model. 

The clustering extends over such large length scales, in 
specific ranges of r/p and t/p^ that the possibility of a 
phase transition into a macroscopic clustered stationary 
state has been the major issue [M, kM]. Macroscopic 
cluster condensation (infinite sized clusters) have been 
shown to be impossible using the analytic matrix product 
ansatz [13 and also using an approximate mapping onto 
the so-called zero range process 0. The cluster size 
remains always finite, but the maximum value can be far 
beyond all computation capabilities [13 • 




FIG. 2: Phase diagram for the AHR model as function of r/p 
and (conserved) average density p — pj^ — p-. 



III. THE PHASE DIAGRAM 

Fig |2] shows the phase diagram of the t = AHR model 
as function of r/p and (conserved) global average density 
p = pj^ = It contains three special lines: r/p = 1, 
r/p = 2^ and p = 0.5, respectively. 

Along the p = 1/2 line all sites are fully occupied and 
the process reduces to the singe-species ASEP. From the 
perspective of the AHR process as modeling two coupled 
conserved degrees of freedom, momentum and density, 
the density sector freezes out, leaving only the Burgers 
equation. The p = 1/2 limit is therefore anomalous, and 
this line is not the proper backbone of the phase diagram. 
The dynamic scaling exponent is equal to z = 3/2 along 
this line, but that does not need to extend to p < 0.5. 

The r/p=l line and the interpretation of the AHR 
process in terms of two coupled Burgers equations form 
the true backbone of the phase diagram. At r = p, the 
process reduces to a single-species ASEP in two different 
ways. If the + particles choose to be blind to the differ- 
ence between an empty site and a — particle, they see 
at r = p no difference between a free hop and a passing 
event, and thus experience pure single species ASEP scal- 
ing. The same is true in the projected subspace where — 
particles are blind to the difference between empty sites 
and + particles. These subspaces are not perpendicular 
and the process does not factorize into two independent 
ASEP processes. Correlations exist between the + and — 
particles, resulting in clustering. We will study this nu- 
merically in the next section and find that at large length 
scales the process factorizes after all, into (KPZ)^. 

At r 7^ p the particles can still pretend to be blind 
to the other species, but then experience updates where 
the hopping probability inexplicably changes from p to 
r. These events are random, but not uncorrelated. For 
r < p the clustering increases and for r > p decreases. 
The line r = 2p is special; there the clustering vanishes 
accidentally altogether. 



IV. DYNAMIC PERFECT SCREENING AT r = p 

Our investigation of dynamic scaling in the AHR model 
started with an attempt to measure the dynamic critical 
exponent z in the conventional manner, i.e., from the 
time evolution of the interface width starting from a fiat 
or a random initial state. Recall that the AHR model is a 
RSOS type KPZ growth model with conserved number of 
up and down steps. It turns out that this interface width 
oscillates in time while evolving toward the stationary 
state, as illustrated in Figl3l 

The fiat initial state evolves roughly in accordance with 
conventional scaling, i.e., as ~ t^, with P = a/z, at 
intermediate times and saturating at ~ (with a 
the stationary state roughness exponent), but the oscil- 
lations on top of this behavior are too strong to accu- 
rately determine /3. These oscillations refiect the addi- 
tional conservation law, and are tied to traveling wave 
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FIG. 3: The evolution of interface widths with uncorrelated 
disordered initial states with r = p = 1.0 and p = 0.25 (a) 
and flat initial states with r — p — 1.^ and p = 1/3 (b) for 
different system sizes L. 



packets propagating in opposite directions and meeting 
again after traveling around the lattice ring. 

For the resolution of this problem we turn our attention 
towards these wave-packets themselves, by monitoring 
the manner they spread in time. This is achieved in terms 
of the two-point correlators 

g+-{x,t) = (n+(0)n_(x))-(n+(0))(n_(x)) (3) 

and and Q--^ where n±{x) is the number operator 
for ± particles at site x and at time t. The perfect screen- 
ing phenomenon in the stationary state emerged while we 
tested this novel method. In this section we first present 
and discuss perfect screening and then present the nu- 
merical analysis of the dynamic exponent, both at r = p. 

A. Stationary State Correlation Functions 

In the stationary state, the correlation function 

g+.{x) = {n+{Q)n.{x)) - p\ (4) 

decays exponentially toward zero. Figs UHS] illustrate this, 
using MC simulations for periodic boundary conditions 
for small rings, L < 800. The correlation function decays 



exponentially for x > and is zero for x < 0. Correla- 
tions are absent at x < 0, because after passing, — and 
+ particles hop away from each other, and (at r = p) do 
not communicate with each other anymore. 

The correlation length is rather short in Fig 21 ^ — 5, 
but increases with density along the r = p line. The most 
significant aspect is not the correlation length, but the 

absence of any finite size scaling offset {x) ^ B/L 

for X ^ ^ and x < 0. 

The absence of this offset is quite surprising. It indi- 
cates a "perfect screening" localization type phenomenon 
in the fluctuations. To appreciate this, consider the two- 
point correlation in a random disordered state, like the 
single species ASEP stationary state. The and Q-- 
correlators in our model have exactly that form at r = p 
because each couples only to one of the two projected 
single species ASEP subspaces. Such correlators are 5- 
functions (with negative 5(0) /L offsets) because periodic 
boundary conditions imply rigorous global conservation 
of the total number of particles, and impose the condi- 
tion that the total area underneath Q{x) is exactly equal 
to zero. 

Another way of viewing this starts by realizing that 

{x)/ p can be interpreted as the probability to find 

a — particle at distance x from a tagged + particle at 
site xq. The tagging removes an amount of probability p 
from xq corresponding to the (untagged) probability of 
finding a — particle at xq. This amount is redistributed 
over the chain. In general, we would expect that part of 
this expelled probability remains localized near xq, rep- 
resented in by the area underneath the exponential; 
and that the remainder is distributed uniformly over the 
chain in delocalized form, represented by a uniform B / L 

type finite size offset in For uncorrelated (^-function 

type correlations all of it is delocalized, such that B = p'^. 
Our numerical simulations, see Fig|5l put a bound on the 
delocalized amplitude; e.g., \B\/L <8x 10"^ at L = 800 
for p = 0.25. The delocalized fraction is zero within the 
MC noise. 

So surprisingly, in our process all the excluded proba- 
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and Q are (^-functions at r = p and their finite size 

offsets are absent in the thermodynamic limit L ^ oo, 

W{x)' = x{p+-pl)+xip.-pl) + 2x{p'-A+_). 

Moreover, at length scales much larger than the screen- 
ing length, X the cross-correlator area = 

"^y^i G-\ {y) reduces to = by perfect screening, 

such that the contributions vanish completely, 

W{xf = x{p+ - pX) ^ x{p_ - pi) 

= W{x',^f ^W{x',-f. (7) 

The square of the full interface width is thus equal to the 
sum of the squared interface widths in the two projected 
subspaces at x^ The two coupled Burgers equations 
behave independently at length scales x ^ 



bility is localized, such that 



C. Dynamic Exponent from {x^t) 



y=l 



(5) 



for all ^ <C a <C i^. A person riding on top of a specific 
+ particle and wearing glasses that filter out the — par- 
ticles, observes a perfect single species ASEP in terms of 
the + particles. Without glasses she notices however an 
excess of — particles in front of her. This cloud of size 
has on average an excess mass equal to p. 



B. Factorization from Perfect Screening 

The above perfect screening implies that the AHR pro- 
cess at r = p behaves at coarse grained scales as two 
decoupled single species Burgers equations. This factor- 
ization is easily recognized in the interface growth repre- 
sentation. Recall that the + particles represent up-steps 
and the — particles down steps, and that the number of 
both are conserved. Perfect screening means factoriza- 
tion into two decoupled KPZ interface growth processes 
at length scales x ^ ^ (one where down-steps are being 
ignored and the other where the up-steps are ignored). 

The interface width W{x) of the full model over a sec- 
tion of the interface of length x can be expressed in terms 
of the two-point correlators as 



W{xf = {{h{x) - h{0)f) 
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Figl6] shows the time evolution of the Q-^ correlation 

function starting from an initial uncorrelated disordered 
state (a (5-function with a finite size off-set). The build- 
up of the cluster of — particles in front of the tagged + 
particle requires only a short time span tq. The build-up 
of this surplus is mirrored by the build-up of a depletion 
layer behind the + particle (particle numbers are locally 
conserved). After the screening cloud at x > is fully 
established, t > tq, the depletion packet detaches from 
X = and travels to the left. This traveling wave packet 
belongs to one of the two projected single species ASEP 
subspaces and therefore should spread in time with KPZ 
dynamic exponent z = 3/2 SiS w ^ t^l^ . The Gaussian 
form 



^a+-(x,t) - t"^/^exp 



{X - Vgif 



(6) 



fits the wave packet very well [2l| except at times close 
to t = To = £,/vg where it is slightly skewed. The 
packet's group velocity, Vg^ follows the expected value 
Vg = 2p{l — 2p), i.e., twice the group velocity of fluc- 
tuations in the — or + sector single species ASEP. {2vg 
is the relative velocity of fluctuations in the + and — 
sectors respectively, propagating in opposite directions.) 
The traveling depletion wave packet moves around the 
ring while broadening. It collides after one period with 
the screening cloud. They split-off again. This keeps re- 
peating itself, until the broadening has spread all over 
the ring and cancels out against the global finite scaling 
offset of the initial state. 

Fig |7l shows the time evolution of the width of the wave 
packet w and its height h. They obey power laws: w ~ 
t^/^ and h ^ t~^l^ . From these, we obtain estimates for 
the dynamic exponent and Fig l8] shows the finite size 
scaling behavior of these estimates. They converge to 
z — 1.53 [2], consistent with the expected KPZ value z — 
3/2. This confirms that this novel way for determining z 
works well. 
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FIG. 6: g+-{x,t) at series of time t =100, 200, 300, 400 with 
r = p = 1.0, p = 0.25, and L = 800. The group velocity is 
equal to Vg — 2p(l — 2p) — 1. 



V. DYNAMIC SCREENING AT r / p 

The correlation functions and 5++ take more in- 
tricate shapes away from the r = p line. Remarkably, as 
we will discuss next, this variety of shapes convert back 
into the simple shapes of r = p using a quasi-particle 
representation. We discovered this numerically, as pre- 
sented in this section, and then proved it analytically, 
as presented in the next section. The properties at the 
r = p line, perfect screening between particles of oppo- 
site charge and uncorrelated disordered stationary state 
statistics in the two projected subspaces, extend thus to 
all r/p in terms of quasi-particles, and the final conclu- 
sion from this is that the process factorizes into (KPZ)^ 
everywhere for all r/p. 
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FIG. 7: Widths and heights of the wave-packets in ^+_(x,t) 
at series of time, t, for p = r = 1.0 and L = 3200. A line 
corresponding to z = 1.5 is drawn above the data points in 
the upper figure. 



A. Stationary State Correlation Functions 

Fig|9] shows the Q-^ and correlators for various 

values of r/p. Compared to the r = p shapes, Q-\ devel- 
ops correlations at x < 0, and = Q changes from 

a (^-function into a symmetric correlated shape. This can 
be explained qualitatively as follows. At r 7^ p, the + 
and — particles can not choose to be bhnd with respect 
to each other anymore. Additional correlations build-up 
compared to the r = p baseline behavior: 

At r < p the passing versus hopping rate is reduced. 

The screening cloud at x > in therefore grows 

(the clustering is stronger). This enhanced Q-\ screen- 
ing cloud at X > 0, results in short range correlations 
between alike particles as well; = Q (x) devel- 
ops positive tails. This is a second order effect. Those 
++ particle correlations in turn induce positive correla- 
tions in Q-i (x) for X < 0. This is a third order effect, 

and thus an order of magnitude further down. 

At r > p the passing rate is enhanced with respect to 
the r = p baseline behavior. The x > screening cloud in 



{x) is thus smaller than at r = p. The correlations in 

= Q are indeed negative, and represent a reduced 

probability to find alike particles near each other. This 
reduced probability makes it less likely to find -h particles 
behind the tagged + particle, at x < 0. If those + parti- 
cles had been there, they would carry smaller screening 
clouds in front of them. Their absence therefore creates 
still positive correlations between — particles at x < 
and the tagged + one. 

At r = 2p the stationary state is fully disordered 
the clustering vanishes and all correlation functions re- 
duce there to (5-functions. At r > 2p the correlation tails 
re-emerge, but with opposite signs. 



B. Dynamic Exponents from {x,t) and ^++(x,t) 

We examine the temporal evolution of Q-\ {x^t) and 

^++(x,t) using MC simulations, just as we did in the 
r = p case. The initial states are prepared to be uncor- 
related and disordered. As shown in FigHOl two wave- 
packets appear, with different amplitudes, but moving 
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in opposite directions with the same speed. The wave- 
packets in Q-\ are strongly coupled to those in 5++. 

These traveling clouds are generated by the same type 
of mechanism as the one at r = i.e., the result of 
the rather fast build-up of the screening clouds near the 
tagged particle, reflected by the short distance correla- 
tions in the stationary state. Both traveling clouds are 
mixtures of + and — particles, with non-zero projections 
in both and G 
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FIG. 10: (a) Correlation function between + and — at 
t = 1000 and 1400 with p = 1.0, r = 0.5, and L = 6400. 
The initial state is random disordered, (b) The correspond- 
ing correlation function between + and +. (c) Q-\ — (x,t) for 
t = 300, 450 with p = 0.7, r = 1.0, and L = 6400. (d) The 
corresponding 5++(x,t). 
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FIG. 11: The estimates of dynamic exponent, z, for different 
system sizes at p = 1.0, r = 0.5, and p = 0.25. 
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FIG. 12: The estimates of dynamic exponent, z, for different 
system sizes ai p = 0.7, r = 1.0, and p = 0.25. 



Once the clouds are detached from x = 0, they move 
independently of each other in opposite directions; just 
as at r = p. The process factorizes again. But there 
is no a priori reason why these mixed traveling clouds 
at r 7^ p should spread as in pure KPZ. However, they 
do. In our MC simulations they spread, e.g., at r = 0.5, 
p 1.0, and p 0.25, with z 1.54[2] and at r = 0.7, 
p = 1.0, and p = 0.25 with z = 1.51[2]. Fig HH and [12] 
shows strong finite size corrections to the scaling in the 
dynamic exponents, but the limiting behavior is clear. 
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FIG. 13: The temporal evolution of the interface widths start- 
ing from initial flat interfaces for p — 0.5, r = 1.0, and 
p = 0.25. (a) The evolution of the widths for different system 
sizes only shows oscillations for i < 200. (b) The dynamic 
exponent is estimated by measuring the slopes of log-log plot 
of the interface width v.s. time. 
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FIG. 14: The estimate of dynamic exponent, z, for different 
system sizes at p = .5, r = 1.0 and p — 0.25. 



Moreover, at r = 2p, the stationary state is totally 
uncorrelated and disordered (and the temporal evolu- 
tions of the correlation functions therefore do not involve 
traveling wave packets). We can apply the conventional 
method to estimate the dynamic exponent. The tem- 
poral evolution of the interface widths (see Fig. (p!3|) and 
(dH), yields z = 1.51[1]. 



C. Perfect Screening in the Quasi-Particle 
Representation 

If indeed the dynamic exponent retains (KPZ)^ type 
value at r 7^ p as suggested by the above numerical re- 
sults, then there might be a quasi-particle description in 
which the process factorizes at large length scales and in 
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which the fluctuations are perfectly screened, just as at 
r = p. We found such a description, flrst numericahy 
as described here, and then rigorously analytically in the 
following section. This implies the process obeys (KPZ)^ 
scaling everywhere. In terms of quasi-particles the dy- 
namic process fully factorizes into two KPZ processes at 
large coarse-grained scales. 

Consider the stationary state correlation functions in 
FiglQlandfTOl the correlation functions decay to 1/L type ''/p 

flnite size scaling off-sets. The area A underneath for 

X > 0, the area B underneath at x > (equal to the 
same for x < 0) and the area C underneath for x < 
obey empirically the relation B / A = C / for all r 7^ p, 
typically with a numerical accuracy, 1 — B'^ /AC = 0.01%. 
(The areas are measured with respect to the offsets.) 

This special balance in the areas relates to a speciflc 
amount of mixing between + and — particles in the 
clouds, and suggests (a much stronger property) the ex- 
istence of a quasi-particle representation. 



f3n^ 



an- 



(8) 



with n± the number operator for + and — particles and 



B 



C_ 
B 



(9) 



in which the correlation functions Qppix) = Grnm{x) and 
Gpm{x), deflned as 



(10) 



with ji = p^m reduce to the same shapes as the parti- 
cle correlators at r/p = 1 (where Gpp{x) is a (5-function 
and Qpm has only one tail and shows perfect screening 
between quasi-particles of opposite charges). 

The mixing ratio R = P/a varies from R = P/a = dit 
r = p (with Up = n+ and Um = n_); to = l3/a=l when 
Up = nm, and to R = P/a = —1 when Up = —nm- FigfT5l 
shows lines of R from our analytic expression in section 
IVIFi Our numerical results are completely consistent 
with this. The mixing strength increases with density p, 
and becomes indeterminate at the line p = 0.5, where all 
sites are fully occupied. At r/p = 2 the stationary state 
is totally disordered, but R does not vanish since + and 
— remain strongly correlated dynamically [22j]. Both a 
and P go to zero and change sign across the r/p = 2 line. 



VI. PERFECT SCREENING AND THE MATRIX 
PRODUCT STATIONARY STATE STRUCTURE 

In this section, we prove analytically the perfect screen- 
ing of the (quasi-particles) pair correlators, using the 
matrix product ansatz (MPA) structure of the station- 
ary state. The proof applies to all r/p, but for clarity 
we split-up the discussion. First we review briefly the 
general properties of MPA stationary states. Next, we 
present the proof at r = p, and flnally generalize it to all 
r/p in terms of quasi-particles. 




FIG. 15: Contour plot of P/a in a parameter space of r/p 
and p (See Ea. (|i9l) ). 



A. MPA Type Stationary States 

Stationary states of stochastic dynamic processes are 
typically very complex with intricate long-range effec- 
tive interactions between the degrees of freedom (when 
writing the stationary state in terms of effective Gibbs- 
Boltzmann factors). The long-range aspect is important; 
ID driven stochastic processes can undergo non-trivial 
phase transitions, while ID equilibrium degrees of free- 
dom with short-range interactions can not. MPA states 
are linked to equilibrium distributions and therefore lack 
long-range correlations. 

MPA stationary states are of the form [l^, [lH, [13] 



Ps{{ri]) = ^Tr[Gr,Gr, 



(11) 



with in our case = +,0,—. This structure resembles 
closely the transfer matrix formulation of partition func- 
tions in one dimensional (ID) equilibrium statistical me- 
chanics. Consider for example a one dimensional Ising 
model, with spin 6* = ±1 degrees of freedom at sites 



^, that interact as 



(12) 



and with a degree of freedom = 0, ±1 on every bond 
i such that bond energy K{ri) can have three distinct 
values. In that case, the Gn are 2x2 transfer matrices 
and the stationary state probability for the yet unrelated 
stochastic dynamic process is the Ising equilibrium par- 
tition function for a given {ri} conflguration. 



^g-i?({T,S})^ 

{S} 



(13) 
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The normahzation factor 



(14) 



is the canonical partition function of the annealed ran- 
dom bond ID Ising model. The stochastic driven non- 
equilibrium dynamics typically imposes constraints on 
the Ti degrees of freedom. In our dynamic process the 
number of each species of particle, r = ±, is conserved 
independently. This is denoted by * in eq. (p!4|) . 

The T variables do not couple to each other directly in 
eq. (p!3|) ; all correlations between r degrees of freedom are 
mediated by the Ising field Sx - The search for a possible 
MPA structure of the stationary state is therefore the 
search for the existence of a representation in which all 
correlations between the original r degrees of freedom are 
carried by a new auxiliary field and expressed as short- 
range interactions between those new degrees of freedom. 
Those auxiliary degrees of freedom can take any form, not 
just Ising spins, because the rank of the G matrices and 
their symmetries can be arbitrary. For example in our 
case, the rank will be infinite, and the auxiliary field can 
be interpreted as (integer valued) interface type degrees 
of freedom, denoted as Tlx — 

The transfer matrix product structure, eq.(pT]). im- 
plies that those auxiliary degrees of freedom interact by 
nearest neighbor interactions only. This is actually un- 
fortunate, because in short-ranged ID equilibrium sys- 
tems, like eq. (p!3|) . spontaneous broken symmetries and 
phase transitions are impossible. Therefore, master equa- 
tions with MPA stationary states have at best dynamic 
phase transitions with trivial scaling properties (associ- 
ated with an abrupt change in the G-representation). For 
example, MPA representations of directed percolation or 
directed Ising type processes can not exist, because both 
are believed to have transitions with complex scaling di- 
mensions. Still, the MPA method has been proved to be 
a powerful tool, its algebraic structure is very elegant, 
and a surprisingly large class of ID stochastic dynamic 
processes have a MPA type stationary state. 

Boundary conditions play an important role. Eq.([T4j) 
is a canonical partition function, where the number of + 
and — particles are each conserved. Consider instead the 
generating function 



Psiin]) 



= Tr[{z+G++Z-G-+Go)^]=Tr[M^] (15) 



with 



M = z+G+ + z-G-+Go. 



(16) 



This would be the grand canonical partition function of, 
e.g., the above annealed random bond ID Ising model 
in case of periodic boundary conditions. z± are the fu- 
gacities of the t± particles. The equivalence between the 
ensembles in the thermodynamic limit is ensured in the 



equilibrium interpretation, where the details of how the 
particle reservoirs couple to the system does not have to 
be addressed. This is different in the interpretation of the 
MPA as the stationary state of a driven stochastic pro- 
cess. Dynamic processes are very sensitive to boundary 
conditions. For example, a process with open boundary 
conditions and reservoirs at the edges conserves the num- 
ber of particles everywhere inside the bulk, and behaves 
very different from the one where the reservoirs couple 
directly to every site. Not surprisingly therefore, the 
MPA method only applies to the stationary state; the 
introduction of the auxiliary field does not address the 
stochastic dynamics, nor the temporal fluctuations in the 
stationary state. For periodic boundary conditions, as in 
our case, the grand canonical partition function eq. (p!5|) 
represents an ensemble of dynamic systems, each with 
periodic boundary condition systems, and fixed values of 
N- and A^+, weighted with respect to each other by the 
fugacity probability factor. In this sense the ensembles 
are equivalent in the thermodynamic limit. In our dis- 
cussion below we use the grand canonical ensemble. The 
correlation functions for x > are evaluated then as 



G-\B) 



-{B\G+\B){B\G-\B)\ (17) 

in the thermodynamic limit, with \B) and {B\ the right 
and left eigenvectors of the largest eigenvalue of the 
operator M defined in eq. (p!6|) . The correlator at x = 
0, (n+n_), poses somewhat of a problem. It can not 
be expressed as simple as this due to the intrinsic off- 
diagonal character of the above G operators. At r = p 
this is not an issue, because (n+n_) = 0. However, that 
will not be true anymore for the quasi-particles at r 7^ p. 



B. Quadratic Algebra 

The first step in identifying whether the stationary 
state of a stochastic process might have a MPA structure, 
is to insert Eq. (pT]) into the master equation. If lucky, the 
condition of stationarity can be expressed as simple alge- 
braic conditions on the Gr transfer matrices. The MPA 
structure of our process has been studied extensively re- 
cently [ii, [H, [ll, [13, [ii, From those studies we 
know that the three Gr must obey the quadratic alge- 
bra: 



rG+G_ 
pGoG_ 



— x_G+ -h x+G_. 
-xoG+ + x+Go, 
-x-Go + xoG_, 



(18) 



with xo and x± arbitrary yet unspecified parameters. 
These conditions apply to the entire phase diagram, for 
all r/p. The next step is to find explicit representations of 
the Gr that satisfy eq.([T8|), using the freedom in choice of 
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the parameters Xi . In general the rank of the Gr does not 
close, but remains infinite. The rank is finite only along 
special lines in the phase diagram. Fortunately, for our 
purposes we do not need closure; the quadratic algebra 
structure itself is sufficient to prove perfect screening. 

Our process is invariant under simultaneous inversion 
in space x —x and of charge + ^ — in the case that 
the numbers of + and — particles are balanced. This sug- 
gests we look for a realization of the algebra with opera- 
tors satisfying G+ = and Go = Gq. This invariance 
is valid in the subspace x+ = — x_ = r and xq = [13], 
where the quadratic algebra reduces to 



Bound states have the generic form 
1 



zw? for n > 0. 



(24) 



Substitution in Eq. (|23|) yields only one bound state, with 
energy Eb = "2 — 1/wb — such that 



= — (1 + ^5) , 



(25) 



and normalization 



G+G- 
G+Go 



G+ + G- 
GqG- = 



Zb 



l-wl 



(26) 



r 

-Go. 

P 



(19) 



Wb is equal to Wb = z. 

The extended eigenstates are scattered waves, with 



C. The r = p Quadratic Algebra 

At r = p, the quadratic algebra of Eq.dTOl) is easily 
checked to be satisfied by the operators [l5| 

G+ = / + L_, G_=/ + L+, Go = |0)(0|. (20) 

The rank of these matrices is infinite. The auxiliary de- 
grees of freedom are (positive only) integer valued "height 
variables" n = 0, 1, 2, • • • . Go is the projection operator 
onto the n = state, and L± are the raising (lowering) 
operators, L±\n) = |n ± 1). 

We need to determine the eigenvalues of the grand 
canonical transfer matrix, eq.(pT 



^o(fc) 



1 



COS Kn 



{kr 



(27) 



The eigenvalue equations at n > 1 yield the energy spec- 
trum, E(k) = 2(1 — cos/c), with < /c < tt, such that 

\{k) = 2^(1 + cos A;), (28) 

and those at n = 0, 1 yield the phase shift Ok^ 

Ao(fe) = cos(^fe) = zcos(^fe -k) = ^) . (29) 



2z - zEk-l 



The normalization factor 

D 



M = zG^^zG-^Go = z{2I^L^^L_)^Go. (21) |^^|2 ^ ^ |^^|2 ^ 1 ^ Z{k) = \Ao\' ^ -D (30) 



This matrix has several interpretations. It is the trans- 
fer matrix for the equilibrium partition function of a one 
dimensional interface in the presence of a substrate (all 
n < are excluded) with a short range attractive poten- 
tial at n = 0; like a substrate. Such an interface layer 
is thin and non-rough. It is also the time evolution of a 
ID random walker (with x playing the role of time and 
n that of space) in half space, n > and an on-site at- 
tractive interaction at site n = 0. Such a random walker 
is localized. The latter can be presented also as localiza- 
tion of a single quantum mechanical particle hopping on a 
semi-infinite chain with a ^-function attractive potential 
at the first site. 



n=l 



is proportional to the rank of the matrices and thus 
strictly speaking infinite; D will drop out in our calcula- 
tions below. The n = component is easily evaluated: 



z sm k 



D z'^ — 2z cos /c + 1 



(31) 



H = 2I-L. 



1 



Go- 



(22) 



with M = 4.z{l- iH). 

This simple Hamiltonian has one single bound state 
and a continuum spectrum of extended states. The cal- 
culation of the eigen-spectrum is elementary and straight 
forward. The eigenstates, |A) = ((/)o, • • • ), satisfy the 
equations: 



D. Perfect Screening at r = p 

Perfect screening implies that 

00 

Y,G+-{^) = -Q+-{0),- (32) 

x=l 

i.e., that the sum over x > of the correlator Eq. ([T7|) . 

CO 

S = J2g+-{x) (33) 



(2- 



■ 2^n 



-(/>!= E(j)Q 



forn > 1. (23) 



is equal to the right hand side of Eq. (|32|) . 

^2 _ 2 I 



(34) 
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with p = zXp/Xb-, using that the bound state is also an 
eigenstate of G+, Gj^\B) = Xp\B) = (1 + Wb)\B). In our 
specific case the density is simply equal to p = z/(l + z), 
but we like to keep the derivation as generic as possible. 

We need to demonstrate that this sum rule is valid in 
the thermodynamic limit, and track carefully any terms 
that scale as system size L. For example, as discussed 
already in detail above, the sum rule is trivially true for 
periodic boundary conditions, but then does not imply 
perfect screening, because any unscreened surplus can 
be spread out over the entire lattice in the form of a 1/L 
background density. 

Define Pb = \B){B\^ as the projection operator onto 
the bound state, and rewrite Eq. (pTj) . as 



x=l 



G-\B). 
(35) 

The bound state does not contribute to the correlators 
inside the sum. Therefore we can project out the bound 
state from M and then perform the summation 



{B\G. 



It 



1 



M 
Xb 



G-\B) 
(36) 



(The single Pb outside the summation originates from 
the X = 1 contribution.) We can remove and G- 
from the above equation, because the bound state is also 
an eigenstate of the lowering operator, G-^\B) = Xp\B)^ 



(0| 



1-^ + Pb 



Pb 



|0) (37) 



writing zG± = M — Gq — zGzp^ and using that Go = |0)(0| 
is the projection operator onto the first site, and also that 
M — XbPb has no projection onto \B). 

The sum rule we seek is now reduced to the identity 



(38) 



withAo = |(B|0)|2 = {B\Go\B) = {B\M-z{G++G-)\B) 
= XB — '^zXp. The left hand side is easily evaluated, using 



Eq. pT]) and that Xb — Xk 



E 



l(o|fe)P 

Xb — Xk 



dk 







— 2zcosk + 1, 

z'^ sin^ k 
(^2 - 22; cos A: + 1)^ * 



(39) 



This is an elementary contour integral along the unit cir- 
cle in the complex w = e*^ plane, with a double pole at 
w = Wb = z within the circle in addition to a single pole 
at = 0. The integral is indeed equal to z'^ / {1 — z^), 
the right hand side of Eq.([37D (A^ = Ab = (1 + zf and 



Ao = l-^'). 



E. Quadratic Algebra at r 7^ p 

The proof of perfect screening for general r/p follows 
the same pattern as at r = p. The operators obeying 
the quadratic algebra conditions, Eq. (fl9|) . are again ex- 
pressed in terms of raising and lowering operators L±, 



G_ 
Go 



+ L_ + (a - l)Go + (s - l)GoL-], 
-[/ + L+ + (a - l)Go + (s - l)L+Go], 



(40) 



where a = r/p and = 1 — (a — 1)^. The transfer matrix 
retains its form 



M = zG+ + zG- + Go = 



with modified Hamiltonian, 



Az 



1 



-H 



(41) 



H 



2/ - L+ - L_ 
-(8-l)(GoL_ 



z 

+ L+Go). 



(42) 



and with 1/z = a/z -\- 2{a — 1). This is again a one- 
dimensional single particle hopping process in a half- 
space, n = 0, 1,2,---. Compared to Eq. ([22|) for a = 
- = 1, the attractive potential at site n = 1 deepens for 
r > p (reducing the clustering and correlation lengths). 
The novel element is the modified hopping probability s 
between sites n = 1 and n = 0. There is still only one 
bound state 



1 



1 1 



for n > 1, 



VZb 



(43) 



with the same functional form for the bound state energy 
as before, Eq. ([25|) . 

The derivation of the extended states is also straight 
forward. They are again of the form, Eq ([27|) . with the 
same energies, Ej. = 2(1 — cos /c), but satisfying the mod- 
ified relations 



sAo{k) = cos{Ok) 



cos{Ok + k) 



2 cos k - 

This leads after some algebra to 

2 _ {w^Wb - l)sin^ A: 2/D 



(wb ^ w^^ — 2 cos k) (w'^ + w\ 



/-I 



2 cos k) 



(44) 



(45) 



v'^ is the second root of the bound state equation Eq.(j33|). 
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F. Quasi-Particle Representation 

We can now identify the exact form of the ratio a/ (5 m 
the quasi-particle representation, eq.®. The represen- 
tation mixes the G± operators in Eq. (|4Q|) as 

Gp = aG+ + (5G- , Gm = + aG_ . (46) 

The projection operator, Go = |0)(0|, and the transfer 
matrix M are invariant. The latter imphes a + /3 = 1. 

The quasi-particle two-point correlation functions take 
the same form as the particle correlators at r = p. In 
particular, the quasi-particle correlation function Qpm{x) 
is zero for all x < 0. This is true when the bound state 
is also an eigenstate of Gp: 

Gp\B) = Xp\B). (47) 

Inserting the bound state, Eq. (|43|) . yields 

aA„ = 1 + awb + — (48) 

and 

^ ^ js' - l)w, -{I -a) 

a {I - a) ' ^ ^ 

The lines of constant (3/ a are shown in Fig. ([15]). (Insert 
the above equations for cc;^, and the relation between a 
and s.) The contours coincide with our numerical results. 

G. Perfect Screening at r ^ p 

The final step is to prove perfect screening in terms of 
the quasi-particles: 

CO 

s = Y,gpm{x) = -gp^{o). (50) 

The left hand side reduces to exactly the same form as 
Eq. ([37|) . using the exact same steps, because the bound 
state is an eigenstate of Gp just like the particle opera- 
tor G+ at r = p; that is all we used there. The right 
hand side is different, because {npUm) = 2aPp is not 
zero anymore. Since a + ^ = 1, it is still true that 
Pp = Pm = P+ = P- = zXp/Xb- Therefore the sum 
rule equation, Eq. ([38|) now takes the form 

with, as before, Aq = K^lO)^ = {B\Go\B) = {B\M - 
z{Gp-\-Gm)\B) = XB — 2zXp. The summation on the left 
leads again to di w = e*^ type contour integral. It still 
has only two poles inside the unit circle: one double pole 
dit w = Wb and one single pole dit w = (with w'^ the 

second root of Eq. (|^3|) .) The result is indeed equal to the 
right hand side after inserting the proper expressions for 
the various eigenvalues and some not very pretty algebra. 



VII. RESULTS AND CONCLUSIONS 

We have studied the two-species asymmetric exclusion 
process (ASEP) to determine whether the addition of a 
local conservation law changes the dynamic scaling prop- 
erties. In the Burgers (hydrodynamics) context the pro- 
cess conserves both momentum and density. In the KPZ 
context it represents interface growth where the numbers 
of up and down steps are conserved. In the ASEP context 
the particle numbers of both species are conserved. 

We find that the dynamic scaling exponent retains the 
KPZ z = 3/2 value. The AHR process factorizes at 
scales larger than the clustering length scale, ^, into two- 
independent KPZ processes. At r = p, where the passing 
and hopping probabilities are equal, this factorization oc- 
curs in terms of + and — particles, while at r 7^ p it is 
established in terms of quasi-particles. This factoriza- 
tion expresses itself as perfect screening between the two 
species of quasi-particles. ^, the screening length, coin- 
cides with the clustering length scale and represents the 
crossover length scale between single KPZ scaling (within 
each cluster) and factorized (KPZ)^ type scaling. 

The conventional method for measuring the dynamic 
exponents in simulations in terms of the time evolution 
of the interface width fails in this process due to the 
presence of time-oscillations with a period proportional 
to the system size; quasi-particles fluctuations have non- 
zero and opposite drift velocities. Instead, we determined 
the dynamic scaling from the two-point correlation func- 
tions. This might be the first time that it is done in this 
manner. 

The stationary state of this process has been studied 
extensively in the recent literature, because it obeys the 
so-called matrix product ansatz (MPA). We used this to 
prove rigorously the factorization of the fluctuations in 
terms of quasi-particles. This previously unknown fea- 
ture of the algebraic structure of the MPA method needs 
to be understood better, in particular in the context of 
clustering phenomena in general. 

The perfect screening phenomenon is clearly a topo- 
logical feature. The above presentation only partially 
exposes those topological properties; by bringing the per- 
fect screening condition into the form of Eqs.(|38|) and 
([5T]) . The right hand sides of both equations only in- 
volve bound state properties. Their left hand sides how- 
ever involve a summation over all extended states; i.e., 
their projections onto n = (|(0|/c)p). The poles of the 
contour integral links this to the bound states and the 
quasi-particle mixing. The formulation of a general proof 
is important because, if topological, the prefect screening 
and (KPZ)^ scaling at large length scales will be generic 
features, valid to many more processes with clustering. 
Its limitations can teach us when and how novel type 
dynamic scaling sets in. 
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